Dynamical simulation of spin-glass and chiral-glass orderings 
in three-dimensional Heisenberg spin glasses 
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Spin-glass and chiral-glass orderings in three-dimensional Heisenberg spin glasses are studied with 
and without randaom magnetic anisotropy by dynamical Monte Carlo simulations. In isotropic 
case, clear evidence of a finite-temperature chiral-glass transition is presented. While the spin 
autocorrelation exhibits only an interrupted aging, the chirality autocorrelation persists to exhibit 
a pronounced aging effect reminisecnt of the one observed in the mean-field model. In anisotropic 
case, asymptotic mixing of the spin and the chirality is observed in the off-equilibrium dynamics. 
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Recently, there arose a growing interest both theoret- 
ically and experimentally in the off-equilibrium dynam- 
ical properties of glassy systems. In particular, aging 
phenomena observed in spin glasses [1] have attracted 
attention of researchers [2]. Unlike systems in thermal 
equihbrium, relaxation of physical quantites depends not 
only on the observation time t but also on the waiting 
time ttu, i.e., how long one waits at a given state before 
the measurements. Recent studies have revealed that the 
ofF-equilibrium dynamics in the spin-glass state generally 
has two characteristic time regimes [2,3]. One is a short- 
time regime, to << t « (to is a microscopic time 
scale), called 'quasi-equilibrium regime', and the other 
is a long-time regime, t >> t^, called 'aging regime' 
or 'out-of-equilibrium regime'. In the quasi-equilibrium 
regime, the relaxation is stationary and the fluctuation- 
dissipation theorem (FDT) holds. The autocorrelation 
function at times and t + 1^ is expected to behave as 
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where is the equilibrium Edwards- Anderson order 
parameter. In the aging regime, the relaxation becomes 
non-stationary, FDT broken, and the autocorrelation 
function decays to zero as i — > oo for fixed i^. 

On theoretical side, both analytical and numerical 
studies of off-equilibrium dynamics of spin glasses have 
so far been limited to Ising-like models, including the 
Edwards- Anderson (EA) model with short-range inter- 
action [4-6] or the mean-field models with long-range in- 
teraction [3,7,8]. Although these analyses on Ising-like 
models succeeded in reproducing some of the features of 
experimental results, many of real spin-glass magnets are 
Heisenberg-like in the sense that the magnetic anisotropy 
is much weaker than the isotropic exchange interaction. 
Thus, in order to make a direct link between theory and 
experiment, it is clearly desirable to study the dynamical 
properties of Heisenberg-like spin-glass models. 

Even at the static level, nature of the experimentally 
observed spin-glass transtition and the spin-glass state is 
not fully understood. Although experiments have pro- 
vided strong evidence that spin-glass magnets exhibit an 



equilibrium phase transition at a finite temperature, nu- 
merical studies indicated that the standard spin-glass 
order occurred only at zero temperature in a three- 
dimensional (3D) Heisenberg spin glass [9-12]. While 
weak magnetic anisotropy inherent to real materials is 
often invoked to explain this apparent discrepancy, it re- 
mains puzzling that no detectable sign of Heisenberg-to- 
Ising crossover has been observed in experiments which is 
usually expected to occur if the observed spin-glass tran- 
sition is caused by the weak magnetic anisotropy [9,10]. 

In order to solve this apparent puzzle, a chirality mech- 
anism of experimentally observed spin-glass transitions 
was recently proposed by the author [11], on the assump- 
tion that an isotropic 3D Heisenberg spin glass exhib- 
ited a finite-temperature chiral-glass transition without 
the conventional spin-glass order, in which only spin- 
reflection symmetry was broken with preserving spin- 
rotation symmetry. 'Chirality' is an Ising-like multi- 
spin variable representing the sense or the handedness 
of the noncoUinear spin structures. It was argued that, 
in real spin-glass magnets, the spin and the chirality 
were "mixed" due to the weak magnetic ansitroy and the 
chiral-glass transition was then "revealed" via anomaly 
in experimentally accessible quantities. Meanwhile, the- 
oretical question whether there really occurs such finite- 
temperature chiral-glass transition in an isotropic 3D 
Heisenberg spin glass, remains somewhat inconclusive 
[11,12]. 

In view of the absence of off-equilibrium simulation 
of Heisenberg spin glasses, and also of the possible im- 
portant role played by the chirality, I will report in the 
present Letter the results of extensive dynamical Monte 
Carlo simulations on isotropic and anisotropic 3D Heisen- 
berg spin glasses, in which the properties of both the spin 
and the chirality are studied. 

The model is the classical Heisenberg model on a sim- 
ple cubic lattice with the nearest-neighbor random Gaus- 
sian couplings, Jij and D'^'^ , defined by the Hamiltonian 
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where Sj ={Sf , Sf , S^) is a three-component unit vec- 
tor, and the sum runs over all nearest-neighbor pairs 
with N = L X L X L spins. Jy is the isotropic ex- 
change coupling with zero mean and variance J, while 
D^^ (/i, p = x,y, z) is the random magnetic anisotropy 
with zero mean and variance D which is assumed to be 
symmetric and traceless, £>,7 = D-f and D^f = 0. 

The local chirality at the i-th site and in the /x-th di- 
rection, Xifii niay be defined for three neighboring spins 
by the scalar [9,12], 
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where (p, 



X, y, z) denotes a unit lattice vector along 
the /i-axis. Note that the chirality defined by Eq.(3) is a 
pseudoscalar in the sense that it is invariant under global 
spin rotation but changes sign under global spin reflec- 
tion or inversion. 

The spin and chirality autocorrelation functions are 
defined by 

Cs{tn„t + t^) = ^Ei[< Si{t^)-Si{t + t^) >], (4) 
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< XiiJ.{tw)XiiJ.{t + tw) >], (5) 



where < • • • > represents the thermal average and [• • •] 
represents the average over bond disorder. 

Monte Carlo simulation is performed based on the 
standard single spin-flip heat-bath method. Starting 
from completely random initial configurations, the sys- 
tem is quenched to a working temperature. Total of 
about 3 X 10^ Monte Carlo steps per spin [MCS] are gen- 
erated in each nm. Sample average is taken over 30-120 
independent bond realizations. The lattice size mainly 
studied is L = 16 with periodic boundary conditions, 
while in some cases lattices with L = 12 and 24 are also 
studied. 

Let us begin with the fully isotropic case, D = 0. The 
spin and chirality autocorrelation functions at a low tem- 
perature T/J = 0.05 are shown in Fig.l as a function of 
t. For larger tyj, the curves of the spin autocorrelation 
function Cs come on top of each other in the long-time 
regime, indicating that the stationary relaxation is recov- 
ered and aging is interrupted. This behavior has been ex- 
pected because the 3D Heisenberg spin glass is believed 
to have no standard spin-glass order [9-12]. Similar in- 
terrupted aging was observed in the 2D Ising spin glass 
which did not have an equilibrium spin-glass order [5] . By 
contrast, the chiral autocorrelation function shows an 
entirely different behavior: FoUwoing the initial decay, it 
exhibits a clear plateau around t ^ and then drops 
sharply for t > t^- It also shows an eminent aging ef- 
fect, namely, as one waits longer, the relaxation becomes 
slower and the plateau-like behavior &t t ^ tw becomes 
more pronounced. 
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(b) chirality 
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FIG.l Spin (a) and chirality (b) autocorrelation func- 
tions of a 3D isotropic Heisenberg spin glass at a tem- 
perature T/J = 0.05 plotted versus log^g^ various 
waiting times t^- The lattice size is i = 16 averaged 
over 66 samples. 
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(b) chirality 
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FIG. 2 The same data as in Figl, but plotted versus 
logio 

In Fig. 2, Cs and are replotted as a function of the 

scaled time t/tw- Reflecting its interrupted aging, the 
curves of Cs for larger t^, now lie below the ones for 
smaller {subaging). By contrast, the curves of for 
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various cross around t/tw ~ 1, and at t > tw, the 
data for larger lie above the ones for smaller tw [su- 
peraging). Such superaging behavior of means that 
the aging in chirality is more enhanced than the one ex- 
pected from the naive t/tw-scaling. Note that, although 
the chirality is an Ising-like variable from symmetry, the 
observed superaging behavior is in contrast to the aging 
behavior of the 3D EA model which was found to satisfy 
a good i/f^u-scaling in the aging regime [5]. It should also 
be noticed that the plateau-like behavior observed here 
has been hardly noticeable in simulations of the 3D EA 
model. Rather, the behavior of C'^ observed here is rem- 
iniscent of the one observed in the mean-field model such 
as the Sherrington-Kirkpatrick (SK) model [7,8]. This 
correspondence might suggest that an effective interac- 
tion between the chiralities is long-ranged. 

While the plateau-like behavior observed in is al- 
ready suggestive of a nonezero chiral Edwards- Anderson 
order parameter, (j^q > 0, more quantitative analysis 
similar to the one recently done by Parisi et al for the 
4D Ising spin glass [6] is performed to extract Qqq from 
the data of in the quasi-equilibrium regime. Finite- 
ness of is also visible in a log-log plot of versus 
t as shown in the inset of Fig. 3, where the data show a 
clear upward curvature. I extract Qqq by fitting the data 
of for i^, = 3 X 10^ to the power-law form of eq.(l) 
in the time range 40 < i < 3,000. The obtained q^Q, 
plotted as a functon of temperature in Fig. 3, clearly in- 
dicates the occurrence of a finite-temperature chiral-glass 
transition at Tcg/J = 0.157 ± 0.01 with the associated 
order-parameter exponent /3cg — 1.1 ± 0.1. The size de- 
pendence turns out to be rather small, although the mean 
values of Qqq tend to slightlty increase around Tea with 
incrcaing L. Since both finite-size effect and fmite-t-uj 
effect tend to underestimate g^Q, one may regard the 
present result as a rather strong evidence of the occur- 
rence of a finite-temperature chiral-glass transition. 
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FIG. 3 Temperature dependence of the Edwards- 
Anderson order parameter of the chirality of a 3D 
isotropic Heisenberg spin glass. The data are averaged 
over 30-120 samples. Inset exhibits the log-log plot of 
the t-dcpcndcncc of the chiral autocorrelation function 
in the quasi-equilibrium regime for L = 16, T/ J = 0.05 
and = 3 X 10^. 



The associated exponent /?cg ^ 1-1 is considerably 
larger than the value of the 3D EA model (3 ^ 0.5 [10], 
but is close to the value of the mean- field model (3 — 1. 
This suggests that the universality class of the chira- 
glass transition of the 3D Heisenberg spin glass might be 
diflferent from that of the; standard 3D Ising spin glass. 
According to the chirality mechanism, the criticality of 
real spin-glass transitions should be the same as that of 
the chiral-glass transition of an isotropic Heisenberg spin 
glass, so long as the magnitude of random anisotropy is 
not too strong. If one tentatively accepts this scenario, 
the present result opens up a new interesting possiblity 
that the universality class of many of real spin-glass tran- 
sitions might differ from that of the standard Ising spin 
glass, contrary to common belief. 

In the presence of weak anisotropy D > 0, chirality 
scenario predicts at the static level that the transition 
behavior of chirality remains essentially the same as in 
the isotropic case, whereas the spin is mixed into the 
chirality, asymptotically showing the same transition be- 
havior as the chirality [11]. In order to see whether such 
"spin-chirality mixing" occurs in the off-equilibrium dy- 
namics, further dynamical simulations are performed for 
the models with random anisotropics D/J = 0.01 ^ 1. 
While chirality exhibits essentially the same dynamical 
behavior as in the isotropic case (not shown here), the be- 
havior of spin at t > t^, changed significantly in the pres- 
ence of anisotropy. As an example, the spin autocorrela- 
tion in the case of weak anisotropy D/ J = 0.01 is shown 
in Fig.4. Even for such small anisotropy, spin is found 
to show superaging behavior asymptotically at t >> 
similar to that of the chirality in the fully isotropic case, 
demonstrating the spin-chirality mixing. 



0.8 



0.2 



spin 


t„=200 « 
t„= 1,000 • 
t„=5,000 n 
. t„=25,000 ■ 
%t t„=1 00,000 ' 


Cs(tw;t+t.) 




D/J=0.01 


V ■ 


T/J=0.05 


L=16 


logio(rtw) 



-4-3-2-10123 

FIG.4 Spin autocorrelation function of the weakly 
anisotropic 3D Heisenberg spin glass with D/J = 0.01 
plottedd versus logiQ{t/tyj). The lattice size is L = 16 
averaged over 60 samples and the temperature is T/ J = 
0.05. 

Experimentally, thcrmoremanent magnetization 
(TRM) or zero-field-cooled (ZFC) magnetization is found 
to show an approximate f/t^-scaling in the aging regime, 
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with small deviation from the perfect scaling in the di- 
rection of subaging [2] . Although this seems in apparent 
contrast to the present result, it should be noticed that 
standard aging experiments have been made by mea- 
suring the magnetic response, not the autocrrelation. 
Recent numerical simulation by Yoshino et al revealed 
that, at least in the case of the SK model, TRM showed 
the subaging even when the spin correlation showed the 
supcraging [13]. Thus, I also calculate the ZFC mag- 
netization for an anisotropic model with D/J = 0.05: 
After the initial quench, the system is evolved in zero 
field during MCS. Then, an external field of instcn- 
sity H / J — 0.05 is turned on and the subsequent growth 
of the magnetization M{t;t^) is recorded. As can be 
seen from Fig. 5, the data show the near t/t^-scaling in 
the aging regime t > where the spin-autocorrelation 
shows the supcraging. Thus, the observed tendency is 
roughly consistent with experiments. It might be inter- 
esting to experimentally investigate the aging properties 
of spin correlations of Heiscnbcrg-like magnets in search 
for possible supcraging behavior. 
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FIG. 5 Zero-field-cooled magnetization of an anisotropic 

3D Heisenberg spin glass with D/ J = 0.05 plotted versus 
logiQt. The field is h/J = 0.05 and the temperature is 
T/ J = 0.05. The lattice size is L = 16 averaged over 80 
samples. 

In summary, equilibrium and off-equilibrium proper- 
ties of the spin and chirality order in 3D Heisenberg spin 
glasses are studied with and without random anisotropy 
by dynamical Monte Carlo simulations. The results are 
basically consistent with the chirality mechanism: In the 



isotropic case, spin and chirality show very different dy- 
namical behaviors consistent with the 'spin-chirality sep- 
aration', whereas in the anisotropic case, spin shows the 
same asymptotic behavior as chirality, consistent with 
the 'spin-chirality mixing' due to magnetic anisotropy. 
Furthermore, clear evidence for the occurrence of a finite- 
temperatmc chiral-glass transition in an isotropic 3D 
Heisenberg spin glass is presented. 
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percomputer center, ISSP, University of Tokyo. 
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